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I. INTRODUCTION 

This contribution is relative to the opening lectures of the ISSAOS 2001 summer school 
and it has the aim to provide the reader with some concepts and techniques concerning 
chaotic dynamics and transport processes in fluids. Our intention is twofold: to give a self- 
consistent introduction to chaos and diffusion, and to offer a guide for the reading of the 
rest of this volume. 

In the following Section we present some basic elements of the chaotic dynamical systems 
theory, as the Lyapunov exponents and the Kolmogorov-Sinai entropy. The third Section is 
devoted to Lagrangian chaos in fluids. The last Section contains an introduction to diffusion 
and transport processes, with particular emphasis on the treatment of non-ideal cases. 



II. SOME BASIC ELEMENTS OF DYNAMICAL SYSTEMS 

A dynamical system may be defined as a deterministic rule for the time evolution of state 
observables. Well known examples are the ordinary differential equations (ODE) in which 
time is continuous: 

^ = f(x(t)), x,fGlR^; (1) 
and maps in which time is discrete: 



x(t + l) = g(x(t)), x,geIR''. (2) 

In the case of maps, the evolution law is straightforward: from x(0) one computes x(l), 
and then x(2) and so on. For ODE's, under rather general assumptions on f, from an initial 
condition x(0) one has a unique trajectory x(t) for t > [47]. Examples of regular behaviors 
(e.g. stable fixed points, limit cycles) are well known, see Figure 1. 
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FIG. 1. Examples of regular attractors: fixed point (above) and limit cycle (below). 

A rather natural question is the possible existence of less regular behaviors i.e. different 
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from stable fixed points, periodic or quasi-periodic motion. 

After the seminal works of Poincare, Lorenz and Henon (to cite only the most eminent 
ones) it is now well established that the so called chaotic behavior is ubiquitous. As a 
relevant system, originated in the geophysical context, we mention the celebrated Lorenz 
model [40]: 

— — —xz + rx — y [6) 
at 

dz , 

— = xy — bz 
dt ^ 

This system is related to the Rayleigh-Benard convection under very crude approximations. 
The quantity x is proportional the circulatory fluid particle velocity; the quantities y and 
z are related to the temperature profile; a, b and r are dimensionless parameters. Lorenz 
studied the case with u = 10 and 6 = 8/3 at varying r (which is proportional to the Rayleigh 
number). It is easy to see by linear analysis that the fixed point (0, 0, 0) is stable for r < 1. 
For r > 1 it becomes unstable and two new fixed points appear: 



C+,_ = {±^b{r - 1), ±^b{r-l), r - 1), (4) 

these are stable for r < Vc — 24.74. A nontrivial behavior, i.e. non periodic, is present for 
r > Tc, as is shown in Figure 2. 

In this "strange", chaotic regime one has the so called sensitive dependence on initial 
conditions. Consider two trajectories, x(i) and x'(t), initially very close and denote with 
A{t) — ||x'(i) — x{t)\\ their separation. Chaotic behavior means that if A(0) — > 0, then as 
t ^ oo one has A{t) ~ A(0) cxp Ait, with Ai > 0, see Figure 3. 

Let us notice that, because of its chaotic behavior and its dissipative nature, i.e. 

dx ^ dy ^ dz ^ ^ . , 

dx dy dz 

the attractor of the Lorenz system cannot be a smooth surface. Indeed the attractor has a 
self-similar structure with a fractal dimension between 2 and 3. The Lorenz model (which 



had an important historical relevance in the development of chaos theory) is now considered 
a paradigmatic example of a chaotic system. 
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FIG. 2. Example of aperiodic signal: the x variable of the Lorenz system (3) as function of time 
t, for r = 28. 
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FIG. 3. The growth of a generic infinitesimal perturbation A{t) in the Lorenz model (3) as 
function of time t, for r = 28. The dashed line represent the exponential growth A(0)expAit. 
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A. Lyapunov exponents 



The sensitive dependence on the initial conditions can be formahzed in order to give it 
a quantitative characterization. The main growth rate of trajectory separation is measured 
by the first (or maximum) Lyapunov exponent, defined as 

Ai = hm hm -In— ^-r, (6) 

t^ooA(0)-»0 t A(0) ^ ^ 

As long as A{t) remains sufficiently small (i.e. infinitesimal, strictly speaking), one can 
regard the separation as a tangent vector z{t) whose time evolution is 

dzi ■^—^ dfi I , , 

-17 = 2^1^Ut)-Zj, (7) 



dt dx,. 



and, therefore. 



A,= limllnMM. (8) 

t^oot ||z(0)|| ^ ^ 

In principle, Ai may depend on the initial condition x(0), but this dependence disappears for 
ergodic systems. In general there exist as many Lyapunov exponents, conventionally written 
in decreasing order Ai > A2 > A3 > as the independent coordinates of the phase space [8]. 
Without entering the details, one can define the sum of the first k Lyapunov exponents as 
the growth rate of an infinitesimal A;— dimensional volume in the phase space. In particular, 
Ai is the growth rate of material lines, Ai + A2 is the growth rate of 2D surfaces, and so on. 
A numerical widely used efficient method is due to Benettin et al. [8]. 

It must be observed that, after a transient, the growth rate of any generic small per- 
turbation (i.e. distance between two initially close trajectories) is measured by the first 
(maximum) Lyapunov exponent Ai, and Ai > means chaos. In such a case, the state of 
the system is unpredictable on long times. Indeed, if we want to predict the state with 
a certain tolerance A then our forecast cannot be pushed over a certain time interval Tp, 
called predictability time, given by: 

Tp ~ ^ In (9) 
Al A(0) ^ ^ 
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The above relation shows that Tp is basically determined by 1/Ai, seen its weak dependence 
on the ratio A/A(0). To be precise one must state that, for a series of reasons, relation (9) 
is too simple to be of actual relevance [14]. 

B. The Kolmogorov-S entropy 

Deterministic chaotic systems, because of their irregular behavior, have many aspects 
in common with stochastic processes. The idea of using stochastic processes to mimic 
chaotic behavior, therefore, is rather natural [21,9]. One of the most relevant and successful 
approaches is symbolic dynamics [7]. For the sake of simplicity let us consider a discrete 
time dynamical system. One can introduce a partition A of the phase space formed by N 
disjoint sets Ai,...,An. Prom any initial condition one has a trajectory 

x(0)^x(l),x(2),...,x(n),... (10) 

dependently on the partition element visited, the trajectory (10), is associated to a symbolic 
sequence 

X(0) ^ il,i2,...,in,- (11) 

where in {in G (1) 2, A'")) means that x(n) e A^^ at the step n, for n — 1,2,.... The 
coarse-grained properties of chaotic trajectories are therefore studied through the discrete 
time process (11). 

An important characterization of symbolic dynamics is given by the Kolmogorov-Sinai 
(K-S) entropy, defined as follows. Let Cn — (ii,«2, ■■■,«n) be a generic "word" of size n and 
P{Cn) its occurrence probability, the quantity 

Hn = SUp[- ^ P{Cn) In P{Cn)] (12) 

A r 

On 

is called block entropy of the n-sequences, and it is computed by taking the largest value 
over all possible partitions. In the limit of infinitely long sequences, the asymptotic entropy 
increment 
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n 



(13) 



is called Kolmogorov- Sinai entropy. The difference — Hn has the intuitive meaning of 
average information gain supplied by the {n + 1)— th symbol, provided that the previous n 
symbols are known. K-S entropy has an important connection with the positive Lyapunov 
exponents of the system [47]: 



In particular, for low-dimensional chaotic systems for which only one Lyapunov exponent is 
positive, one has Hks = Ai. 

We observe that in (12) there is a technical difficulty, i.e. taking the sup over all the 
possible partitions. However, sometimes there exits a special partition, called generating 
partition, for which one finds that Hn coincides with its superior bound. Unfortunately 
the generating partition is often hard to find, even admitting that it exist. Nevertheless, 
given a certain partition, chosen by physical intuition, the statistical properties of the related 
symbol sequences can give information on the dynamical system beneath. For example, if the 
probability of observing a symbol (state) depends only by the knowledge of the immediately 
preceding symbol, the symbohc process is called a Markov chain and all the statistical 
properties are determined by the transition matrix elements Wij giving the probability of 
observing a transition i — > j in one time step. If the memory of the system extends far 
beyond the time step between two consecutive symbols, and the occurrence probability of a 
symbol depends on k preceding steps, the process is called Markov process of order k and, 
in principle, a k rank tensor would be required to describe the dynamical system with good 
accuracy. It is possible to demonstrate that if i?n+i — Hn — hxs ior n > k + 1, k is the 
(minimum) order of the required Markov process [32]. It has to be pointed out, however, 
that to know the order of the suitable Markov process we need is of no practical utility if 
A;> 1. 

For applications of the Markovian approach to geophysical systems see [19] and the 
contributions by Abel et al. and by Pasmanter et al. in this volume. 
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Ai>0 



III. LAGRANGIAN CHAOS 



A problem of great interest concerns the study of the spatial and temporal structure of 
the so-called passive fields, indicating by this term passively quantities driven by the flow, 
such as the temperature under certain conditions [45] . The equation for the evolution of a 
passive scalar field ^(x, t), advected by a velocity field v(x, t), is 

dtO + V ■ (ve) ^X^^O (15) 

where v(x, t) is a given velocity field and x is the molecular diffusion coefficient. 

The problem (15) can be studied through two different approaches. Either one deals at 
any time with the field 9 in the space domain covered by the fluid, or one deals with the 
trajectory of each fluid particle. The two approaches are usually designed as "Eulerian" and 
"Lagrangian" , although both of them are due to Euler [34]. The two points of view are in 
principle equivalent. 

The motion of a fluid particle is determined by the differential equation 

f^v(x.O (16) 

which also describes the motion of test particles, for example a powder embedded in the 
fluid, provided that the particles are neutral and small enough not to perturb the velocity 
field, although large enough not to perform a Brownian motion. Particles of this type 
are commonly used for flow visualization in fluid mechanics experiments, see [51] and the 
contribution of Cenedese et al. in this volume. Let us note that the true equation for the 
motion of a material particle in a fluid can be rather complicated [43,22]. 

It is now well established that even in regular velocity field the motion of fluid par- 
ticles can be very irregular [29,2]. In this case initially nearby trajectories diverge expo- 
nentially and one speaks of Lagrangian chaos. In general, chaotic behaviors can arise in 
two-dimensional flow only for time dependent velocity fields in two dimensions, while it can 
be present even for stationary velocity fields in three dimensions. 
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If X = 0, it is easy to realize that (15) is equivalent to (16). In fact, we can write 

e(x,i) = e„(r-*x) (17) 

where ^o(x) — 9{x.,t — 0) and T is the formal evolution operator of (16) , 

x(i)=T*x(0). (18) 

Taking into account the molecular diffusion x, (15) is the Fokker-Planck equation of the 
Langevin equation [20] 

^ = v(x,t)+r;(t) (19) 

where 77 is a Gaussian process with zero mean and variance 

{rj,{t)rjj{t')) = 2xSij5{t-t'). (20) 

In the following we will consider only incompressible flow 

V-v = (21) 

for which the dynamical system (16) is conservative. In two dimensions, the constraint (21) 
is automatically satisfied assuming 

where ip{'x, t) is the stream function. Inserting (22) into (16) the evolution equations become 

dxi _ di^ dx2 _ _ 
dt 8x2^ dt dxi' 

Formally (23) is a Hamiltonian system with the Hamiltonian given by the stream function 



A. Examples of Lagrangian chaos 



As a first example we consider a 3d stationary velocity field, the so-called ABC fiow 

V = {Asinz + C cos y, Bsinx + A cos z, Csiny + B cos x) (24) 

where A, B and C are non zero real parameters. Because of the incompressibility condition, 
the evolution x(0) — > x(t) defines a volume preserving, dynamics. 

Arnold [4] argued that (24) is a good candidate for chaotic motion. Let us briefly repeat 
his elegant argument. For a steady state solution of the 3d Euler equation one has: 

V- v = 

V X (V X v) = Va (25) 

P v2 

a = — + — 
P 2 

where P is the pressure and p the density. As a consequence of the Bernoulli theorem [35] , 
q;(x) is constant along a streamline - that is a Lagrangian trajectory x(t). One can easily 
verify that chaotic motion can appear only if Q!(x) is constant (i.e. Vq!(x) = 0) in a part 
of the space. Otherwise the trajectory would be confined on a 2d surface q;(x) = constant, 
where the motion must be regular as a consequence of general arguments [47]. In order to 
satisfy such a constraint, from (25) one has the Beltrami condition: 

V X V = 7(x) V. (26) 

The reader can easily verify that the field v given by (24) satisfy (26) (in this case 7(x) = 
constant). Indeed, numerical experiments by Henon [29] provided evidence that Lagrangian 
motion under velocity (24) is chaotic for typical values of the parameters A, B, and C (see 
an example in Figure 4). 
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FIG. 4. Intersections with the Poincare section, plane ^; = 0, of eight trajectories of the ABC 
flow with parameters A = 2.0, B = 1.70, C = 1.50. 

In a two-dimensional incompressible stationary flows the motion of fluid particles is given 
by a time independent Hamiltonian system with one degree of freedom and, since trajectories 
follow iso-V' lines, it is impossible to have chaos. However, for explicit time dependent stream 
function ijj the system (24) can exhibit chaotic motion [47]. 

In the particular case of time periodic velocity fields, v(x, i + T) = v(x, t), the trajectory 
of (16) can be studied in terms of discrete dynamical systems: the position x(t + T) is 
determined in terms of x(t). The map x(t) — > x(t + T) will not depend on t thus (16) can 
be written in the form 



where now the time is measured in units of the period T. Because of incompressibility, the 
map (27) is conservative: 



x(n+l) =F[x(n)], 



(27) 



det^[x]| = 1, 



where Aij [x] 




(28) 
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An explicit deduction of the form of F for a general 2d or 3d flow is usually very difficult. 
However, in some simple model of can be deduced on the basis of physical features [3,18]. 



B. Eulerian properties and Lagrangian chaos 



In principle, the evolution of the velocity field v is described by partial differential equa- 
tions, e.g. Navier-Stokes or Boussinesq equations. However, often in weakly turbulent 
situations, a good approximation of the flow can be obtained by using a Galerkin approach, 
and reducing the Eulerian problem to a (small) system of F ordinary differential equations 
[15,37]. The motion of a fluid particle is then described by the ((i+F)-dimensional dynamical 
system 



where d is the space dimensionality and Q = {Qi, ■■■Qf) are the F variables, usually normal 
modes, which are a representation of the velocity fleld v. Note that the Eulerian equations 
(29) do not depend on the Lagrangian part (30) and can be solved independently. 

In order to characterize the degree of chaos, three different Lyapunov exponents can be 
deflned [25]: 

• a) Ae for the Eulerian part (29); 

• b) Al for the Lagrangian part (30), where the evolution of the velocity fleld is assumed 
to be known; 

• c) At per for the total system of the d + F equations. 
These Lyapunov exponents are defined as: 



dQ 



f (Q, t) with Q, f e 



(29) 



Itt 



v(x, Q) with X, V e H'' 



(30) 




(31) 
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where the evolution of the three tangent vectors z are given by the hnearized stabihty equa- 
tions for the Eulerian part, for the Lagrangian part and for the total system, respectively: 



dt p{ dQj 



z/^\ z(E) e (32) 



Q(t) 



dt dxj 

3=1 J 



Z. 



(L), z(L) elR'^ (33) 



c(t) 



E 



zP, z(T) e IR^+'^ (34) 

y(*) 



dt p{ dyj 

and y = {Qi, . . . , Qf, Xi, . . . , x^) and G = (/i, . . . , fp, Vi, . . . ,Vij). The meaning of these 
Lyapunov exponents is evident: 

• a) Ae is the mean exponential rate of the increasing of the uncertainty in the knowledge 
of the velocity field (which is, by definition, independent on the Lagrangian motion); 

• b) Al estimates the rate at which the distance 5x{t) between two fluid particles initially 
close increases with time, when the velocity field is given, i.e. a particle pair in the 
same Eulerian realization; 

• c) At is the rate of growth of the distance between initially close particle pairs, when 
the velocity field is not known with infinite precision. 

There is no general relation between Ae and Al. One could expect that in presence of a 
chaotic velocity field the particle motion has to be chaotic. However, the inequality Al > Ae 
- even if generic - sometimes does not hold, e.g. in some systems like the Lorenz model [25] 
and in generic 2d fiows when the Lagrangian motion happens around well defined vortex 
structures [6] as discussed in the following. On the contrary, one has 



At = max(AE, Al). (35) 
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C. Lagrangian chaos in two dimensional flows 



Let us now consider the two-dimensional Navier-Stokes equations with periodic boundary 
conditions at low Reynolds numbers, for which we can expand the stream function ip in 
Fourier series and takes into account only the first F terms [15,37], 

F 

^^-iY^ kJ^Q,(^^^- + c.c., (36) 

where c.c. indicates the complex conjugate term and Q = (Qi, . . . , Qp) are the Fourier 
coefficients. Inserting (36) into the Navier-Stokes equations and by an appropriate time 
rescaling, we obtain the system of F ordinary differential equations 

^ = -fcj Q, + ^ A.iMQm + /i, (37) 

^'^ l,m 

in which fj represents an external forcing. 

Franceschini and coworkers have studied this truncated model with F = 5 and F = 7 
[15,37]. The forcing were restricted to the 3*'* mode fj = ReSj^s [37]. For F = 5 and 
Re < Rci = 22.85 . . ., there are four stable stationary solutions, say Q, and Ae < 0. At 
Re = Rei, these solutions become unstable, via a Hopf bifurcation [42], and four stable 
periodic orbits appear, still implying Ae = 0. For Rei < Re < Re2 = 28.41 . . ., one thus 
finds the stable hmit cycles: 

Q(t) = Q + (Re - Rei)^/^(5Q(i) + 0(Re - Rei) (38) 

where dQ,{t) is periodic with period 

T(Re) =To + 0(Re-Rei) Tq = 0.7328 .. . (39) 

At Re = Re2, these limit cycles lose stability and there is a period doubling cascade toward 
Eulerian chaos. 

Let us now discuss the Lagrangian behavior of a fluid particle. For Re < Rei, the 
stream function is asymptotically stationary, V'(x, i) — > V'(x), and the corresponding one- 
dimensional Hamiltonian is time-indcpcndcnt, therefore Lagrangian trajectories are regular. 
For Re = Rei -|- e the stream function becomes time dependent 
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V'(x,t) =V'(x) + v^5V'(x,t) + 0(e), (40) 

where ■0(x) is given by Q and Sil^ is periodic in x and in t with period T. The region of 
phase space, here the real two-dimensional space, adjacent to a separatrix is very sensitive to 
perturbations, even of very weak intensity. Figure 5 shows the structure of the separatrices, 
i.e. the orbits of infinite periods at Re = Rci — 0.05. 




FIG. 5. Structure of the separatrices in the 5-mode model (36) with Re = Rei — 0.05. 

Indeed, generically in one- dimensional Hamiltonian systems, a periodic perturbation 
gives origin to stochastic layers around the separatrices where the motion is chaotic, as 
consequence of unfolding and crossing of the stable and unstable manifolds in domains 
centered at the hyperbolic fixed points [21,47]. One has strong numerical evidence for the 
existence of the chaotic regions, see Figure 6. 

Chaotic and regular motion for small e = Rci — Re can be studied by the Poincare map 

x(nr) ^ x(nr + r). (41) 

The period T(e) is computed numerically. The size of the stochastic layers rapidly increase 
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with e. At e = ec ~ 0.7 they overlap and it is practically impossible to distinguish between 
regular and chaotic zones. At e > ec there is always diffusive motion. 




1 2 3 4 5 6 

FIG. 6. Poincare map for three trajectories of the 5-mode model with Re = Rei + 0.05. The 
initial conditions are selected close to a separatrix, case a) (.xi(0) = 3.2, X2(0) = —1.6), or far from 
the separatrices, cases b) (xi(0) = 4.3, 0:2(0) = -2.0) and c) (xi(0) = 4.267, 0:2(0) = -3.009). 

We stress that this scenario for the onset of Lagrangian chaos in two-dimensional fluids is 
generic and does not depend on the particular truncated model. In fact, it is only related to 
the appearance of stochastic layers under the effects of small time-dependent perturbations 
in one-dimensional integrable Hamiltonian systems. As consequence of a general features of 
one-dimensional Hamiltonian systems we expect that a stationary stream function becomes 
time periodic through a Hopf bifurcation as occurs for all known truncated models of Navier- 
Stokes equations. 

We have seen that there is no simple relation between Eulerian and Lagrangian behaviors. 
In the following, we shall discuss two important points: 

• (i) what are the effects on the Lagrangian chaos of the transition to Eulerian chaos. 
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i.e. from Ae = to Ae > 0. 

• (ii) whether a chaotic velocity field (Ae > 0) always implies an erratic motion of fiuid 
particles. 

The first point can be studied again within the F = 5 modes model (37). Increasing Re, 
the limit cycles bifurcate to new double period orbits followed by a period doubling transition 
to chaos and a strange attractor appears at Rcc ~ 28.73, where Ae becomes positive. These 
transitions have no signature on Lagrangian behavior, as it is shown in Figure 7, i.e. the 
onset of Eulerian chaos has no influence on Lagrangian properties. 
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FIG. 7. Lyapunov exponents Ae (+) and Xl (x) as function of Re around Rec, for the 5-mode 
model. 

This feature should be valid in most situations, since it is natural to expect that in generic 
cases there is a strong separation of the characteristic times for Eulerian and Lagrangian 
behaviors. 

The second point - the conjecture that a chaotic velocity field always implies chaotic 
motion of particles - looks very reasonable. Indeed, it appears to hold in many systems [23]. 
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Nevertheless, one can find a class of systems where it is false, e.g. the equations (29), (30) 
may exhibit Eulerian chaoticity Ae > 0, even if Al = [6]. 

Consider for example the motion of N point vortices in the plane with circulations Fj 
and positions {xi{t) , yi{t)) {i — 1, ..N) [1]: 

„ dxi OH , ^. 

dt dyi 

where 

// = -i-^r,r,lnr,, (44) 

and = [xi - Xjf + {yi - yjf. 

The motion of N point vortices is described in an Eulerian phase space with 2N dimen- 
sions. Because of the presence of global conserved quantities, a system of three vortices is 
integrable and there is no exponential divergence of nearby trajectories in phase space. For 

> 4, apart from non generic initial conditions and/or values of the parameters Fj, the 
system is chaotic [1] . 

The motion of a passively advected particle located in {x{t),y{t)) in the velocity field 
defined by (42-43) is given 

dx ^Tj y-yj 

-dt^-^2^^ ^^^^ 

dy \ ^ Fj X Xi f An\ 

-dt=^2^^ ^''^ 

where = {x — XiY + {y — yiY ■ 

Let us first consider the motion of advected particles in a three- vortices (integrable) 
system in which Ae = 0. In this case, the stream function for the advected particle is periodic 
in time and the expectation is that the advected particles may display chaotic behavior. 
The typical trajectories of passive particles which have initially been placed respectively in 
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close proximity of a vortex center or in the background field between the vortices display 
a very different behavior. The particle seeded close to the vortex center displays a regular 
motion around the vortex and thus \l = 0; by contrast, the particle in the background field 
undergoes an irregular and aperiodic trajectory, and Al is positive. 

We now discuss a case where the Eulerian fiow is chaotic i.e. N — A point vortices. 
Let us consider again the trajectory of a passive particle deployed in proximity of a vortex 
center. As before, the particle rotates around the moving vortex. The vortex motion is 
chaotic; consequently, the particle position is unpredictable on large times as is the vortex 
position. Nevertheless, the Lagrangian Lyapunov exponent for this trajectory is zero (i.e. 
two initially close particles around the vortex remain close) , even if the Eulerian Lyapunov 
exponent is positive, see Figure 8. 
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FIG. 8. Particle trajectories in the four-vortex system. Eulerian dynamics in this case is chaotic. 
The left panel shows a regular Lagrangian trajectory while the right panel shows a chaotic La- 
grangian trajectory. The different behavior of the two particles is due to different initial conditions. 

This result indicates once more that there is no strict link between Eulerian and La- 
grangian chaoticity. 

One may wonder whether a much more complex Eulerian fiow, such as 2d turbulence, 
may give the same scenario for particle advection: i.e. regular trajectories close to the 
vortices and chaotic behavior between the vortices. It has been shown that this is indeed 
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the case [6] and that the chaotic nature of the trajectories of advected particles is not strictly 
determined by the complex time evolution of the turbulent flow. 

We have seen that there is no general relation between Lagrangian and Eulerian chaos. In 
the typical situation Lagrangian chaos may appear also for regular velocity fields. However, 
it is also possible to have the opposite situation, with Al in presence of Eulerian chaos, as 
in the example of Lagrangian motion inside vortex structures. As an important consequence 
of this discussion we remark that it is not possible to separate Lagrangian and Eulerian 
properties in a measured trajectory, e.g. a buoy in the oceanic currents [46]. Indeed, using 
the standard methods for data analysis [27], from Lagrangian trajectories one extracts the 
total Lyapunov exponent At and not Al or Ae. 

IV. TRANSPORT AND DIFFUSION 

The simplest model of diffusion is the Brownian motion, the erratic movement of a 
grains suspended in liquid observed by the botanist Robert Brown as early as in 1827. After 
the fundamental work of Einstein [24] and Langevin [36], Brownian motion become the 
prototypical example of stochastic process. 

A. The random walk model for Brownian motion 

In order to study more in detail the properties of diffusion, let us introduce the simplest 
model of Brownian motion, i.e. the one- dimensional random walk. The walker moves on a 
line making discrete jumps = ±1 at discrete times. The position of the walker, started at 
the origin eit t — 0, will be 

R{t) = i2v, (47) 

i=l 

Assuming equiprobability for left and right jumps (no mean motion), the probability that 
at time t the walker is in position x will be 
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(48) 



For large t and x (i.e. after many microscopic steps) we can use Stirling approximation and 

get 



Pt{x) 



t 



27r(i2 - x^) 



exp 



t + x t + x t — X t — X 
In in 



(49) 



2 2 2 2 

The core of the distribution recovers the well known Gaussian form, i.e. lor x <^t from (49) 
we get 

,.2\ 



(50) 



From (50) one obtains that the variance of the displacement follows diffusive behavior, 



I.e. 



{x'{t))=2t 



(51) 



We stress that diffusion is obtained only asymptotically (i.e. for t — > oo). This is a conse- 
quence of central limit theorem which assures Gaussian distributions and diffusive behavior 
in the limit of many independent jumps. The necessary, and sufficient, condition for observ- 
ing diffusive regime is the existence of a finite correlation time (here represented by discrete 
time between jumps) for the microscopic dynamics. Let us stress that this is the important 
ingredient for diffusion, and not a stochastic microscopic dynamics. We will see below that 
diffusion can arise even in completely deterministic systems. 

Another important remark is that Gaussian distribution (50) is intrinsic of diffusion 
process, independent on the distribution of microscopic jumps: indeed only the first two 
moments of Vi enter into expression (50). This is, of course, the essence of the central limit 
theorem. Following the above derivation, it is clear that this is true only in the core of the 
distribution. The far tails keep memory of the microscopic process and are, in general, not 
Gaussian. As an example, in Figure 9 we plot the pdf pt{x) at step t = 100 compared with 
the Gaussian approximation. Deviations are evident in the tails. 
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FIG. 9. Probability distribution function of the one-dimensional random walk after t 
steps. The dashed line is the Gaussian distribution 
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The Gaussian distribution (50) can be obtained as the solution of the diffusion equation 
which governs the evolution of the probability in time. This is the Fokker-Planck equation 
for the particular stochastic process. A direct way to relate the one-dimensional random 
walk to the diffusion equation is obtained by introducing the master equation, i.e. the time 
evolution of the probability [26] : 

(52) 



Pt+i{x) = ^Pt{x - 1) + ^Pt{x + 1) . 



In order to get a continuous limit, we introduce explicitly the steps Ax and At and write 

Pt+At{x) -Ptjx) ^ {AxYpt{x + Ax) +ptix- Ax) - 2pt{x) 

At ~ 2At {Axf ' ^ ' 

Now, taking the hmit Ax, At ^ in such a way that {AxY / At — > 2D (the factor 2 is purely 
conventional) we obtain the diffusion equation 



dp{x, t) ^ ^fp{x, t) 



dt 



dx'^ 



(54) 



The way the limit Ax, At — » is taken reflects the scaling invariance property of diffusion 
equation. The solution to (54) is readily obtained as 
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Diffusion equation (54) is liere written for tlie probability p{x, t) of observing a marked 
particle (a tracer) in position x at time t. The same equation can have another interpre- 
tation, in which p{x, t) — 9{x, t) represents the concentration of a scalar quantity (marked 
fluid, temperature, pollutant) as function of time. The only difference is, of course, in the 
normalization. 



B. Less simple transport processes 

As already stated time decorrelation is the key ingredient for diffusion. In the random 
walker model it is a consequence of randomness: the steps Vi are random uncorrelated 
variables and this assures the applicability of central limit theorem. But we can have a 
finite time correlation and thus diffusion also without randomness. To be more specific, let 
us consider the following deterministic model (standard map [21]): 

J(t + 1) = J(t) + Ksme(t) 

(56) 

e{t + i) = e{t) + j{t + i) 
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FIG. 10. Square dispersion {J{t)^) for the standard map at K = 10.5. The dashed hne is the 
RPA prediction. 

The map is known to display large-scale chaotic behavior for K > Kc ^ 0.9716 and, as 
a consequence of deterministic chaos, J{t) has diffusive behavior. For large times, J{t) is 
large and thus the angle 9{t) rotates rapidly. In this limit, we can assume that at each step 
9{t) decorrelates and thus write 

J{t)^^K''(j2^sm9{t')j '^K''{siv?e)t^2Dt (57) 

The diffusion coefficient in the random phase approximation, i.e. assuming that sin6{t) is 
not correlated with sinO{t') for t ^ t', is obtained by the above expression as Drpa = 
In Figure 10 we plot a numerical simulation obtained with the standard map. Diffusive 
behavior is clearly visible at long time. 

The two examples discussed above are in completely different classes: stochastic for the 
random walk (47) and deterministic for the standard map (56). Despite this difference 
in the microscopic dynamics, both lead to a macroscopic diffusion equation and Gaussian 
distribution. This demonstrates how diffusion equation is of general applicability. 

C. Advection— diffusion 

Let us now consider the more complex situation of dispersion in a non-steady fluid 
with velocity field v(x, t). For simplicity will we consider incompressible flow (i.e. for 
which V • V = 0) which can be laminar or turbulent, solution of Navier-Stokes equation or 
synthetically generated according to a given algorithm. In presence of v(x, t) , the diffusion 
equation (54) becomes the advection- diffusion equation for the concentration 9{'x,t) (15). 
This equation is linear in 9 but nevertheless it can display very interesting and non trivial 
properties even in presence of simple velocity fields, as a consequence of Lagrangian chaos. 
In the following we will consider a very simple example of diffusion in presence of an array 
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of vortices. The example will illustrate in a nice way the basic mechanisms and effects of 
interaction between deterministic (v) and stochastic (D) components. 

Let us remark that we will not consider here the problem of transport in turbulent 
velocity field. This is a very classical problem, with obvious and important applications, 
which has recently attracted a renewal theoretical interest as a model for understanding the 
basic properties of turbulence [48] . 

Before going into the example, let us make some general consideration. We have seen 
that in physical systems the molecular diffusivity is typically very small. Thus in (15) the 
advection term dominates over diffusion. This is quantified by the Peclet number, which is 
the ratio of the typical value of the advection term to the diffusive term 



where Vq is the typical velocity at the typical scale of the flow Iq. With tq ~ lo/t^o we will 
denote the typical correlation time of the velocity. 

The central point in the following discussion is the concept of eddy diffusivity. The idea is 
rather simple and dates back to the classical work of Taylor [50]. To illustrate this concept, 
let us consider a Lagrangian description of dispersion in which the trajectory of a tracer x(i) 
is given by (16). Being interested in the limit Pe — > oo, in the following we will neglect, just 
for simplicity, the molecular diffusivity D, which is generally much lesser that the effective 
dynamical diffusion coefficient. 

Starting from the origin, x(0) = 0, and assuming (v) = we have (x(t)) = for ever. 
The square displacement, on the other hand, grows according to 



where we have introduced, for simplicity of notation, the Lagrangian velocity vz,(t) — 
v{x{t),t). Define the Lagrangian correlation time tl from 




(-x(t)2) = (x(t) • Mt)) = (v^.) • Mt))ds 



(59) 




(60) 
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and assume that the integral converge so that tl is finite. Prom (59), for t tl we get 



(61) 



i.e. diffusive behavior with diffusion coefficient (eddy diffusivity) — tl(v|). 

This simple derivation shows, once more, that diffusion has to be expected in general in 
presence of a finite correlation time tl- Coming back to the advection-diffusion equation (15), 
the above argument means that for i ^ tl we expect that the evolution of the concentration, 
for scales larger than lo, can be described by an effective diffusion equation, i.e. 



dt 



(62) 



The computation of the eddy diffusivity for a given Eulerian fiow is not an easy task. It can 
be done explicitly only in the case of simple flows, for example by means of homogenization 
theory [41,10]. In the general case it is relatively simple [10] to give some bounds, the 
simplest one being > D, i.e. the presence of a (incompressible) velocity field enhances 
large-scale transport. To be more specific, let us now consider the example of transport in 
a one-dimensional array of vortices (cellular fiow) sketched in Figure 11. This simple two- 
dimensional fiow is useful for illustrating the transport across barrier. Moreover, it naturally 
arises in several fiuid dynamics contexts, such as, for example, convective patterns [49]. 



FIG. 11. Cellular flow model. Iq is the size of vortices, 6 is the thickness of the boundary layer. 



Let us denote by vq the typical velocity inside the cell of size and let D the molecular 
diffusivity. Because of the cellular structure, particles inside a vortex can exit only as a 
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consequence of molecular diffusion. In a characteristic vortex time tq ~ Iq/vq, only the 
particles in the boundary layer of thickness 5 can cross the separatrix where 

5^ ~ Dto ~ D— . (63) 

These particles are baUistically advected by the velocity field across the vortex so they see 
a "diffusion coefficient" 11/tq. Taking into account that this fraction of particles is 5/Iq we 
obtain an estimation for the effective diffusion coefficient as 

^ ^ JdIoVq ~ L>PeV2 (64) 

The above result, which can be made more rigorous, was confirmed by nice experiments 
made by Solomon and GoUub [49]. Because, as already stressed above, typically Pe ^ 1, 
one has from (64) that ^ D. On the other hand, this result do not mean that molecular 
diffusion D plays no role in the dispersion process. Indeed, if D = there is not mechanism 
for the particles to exit from vortices. 

Diffusion equation (62) is the typical long-time behavior in generic flow. There exist also 
the possibility of the so-called anomalous diffusion, i.e. when the spreading of particle do 
not grow linearly in time, but with a power law 

(^'(O) ~ t'" (65) 

with V ^ 1/2. The case v > 1/2 (formally — oo) is called super-diffusion; sub-diffusion, 
i.e. u < 1/2 (formally = 0), is possible only for compressible velocity fields. 

Super-diffusion arises when the Taylor argument for deriving (61) fails and formally 
— > oo. This can be due to one of the following mechanisms: 

a) the divergence of (v|) (which is the case of Levy flights), or 

b) the lack of decorrelation and thus ^ oo {Levy walks). The second case is more physical 
and it is related to the existence of strong correlations in the dynamics, even at large times 
and scales. 
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One of the simplest examples of Levy walks is the dispersion in a quenched random shear 
flow [16,31]. The flow, sketched in Figure 12, is a super-position of strips of size 5 of constant 
velocity vq with random directions. 



FIG. 12. Random shear of ibuo velocity in strips of size 5 

Let us now consider a particle which moves according to the flow of Figure 12. Because 
the velocity field is in the x direction only, in a time t the typical displacement in the y 
direction is due to molecular diffusion only 



5y - \/Dt (66) 

and thus in this time the walker visits N = 6y/6 strips. Because of the random distribution 
of the velocity in the strips, the mean velocity in the N strips is zero, but we may expect 
about y/N unbalanced strips (say in the right direction) . The fraction of time t spent in the 
unbalanced strips is t\/N /N and thus we expect a displacement 

5x ~ Vn^= . (67) 

Wn ^ ' 

From (66) we have ~ \fDi/5 and finally 



i.e. a super-diffusive behavior with exponent v — 3/4. 

The origin of the anomalous behavior in the above example is in the quenched nature 
of the shear and in the presence of large stripes with positive (or negative) velocity in the 
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X direction. This leads to an infinite decorrelation time for Lagrangian tracers and thus to 
a singularity in (61). We conclude this example by observing that for D (68) gives 
{Sx^) — > oo. This is not a surprise because in this case the motion is ballistic and the correct 
exponent becomes v — 1. 

As it was in the case of standard diffusion, also in the case of anomalous diffusion the 
key ingredient is not randomness. Again, the standard map model (56) is known to show 
anomalous behavior for particular values of K [52] . An example is plotted in Figure 13 for 
K = 6.9115 in which one find {J{tf) - t^'^^. 
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FIG. 13. Square dispersion {J{tf) for the standard map at if = 6.9155. The dashed line is t^-^^. 

The qualitative mechanism for anomalous dispersion in the standard map can be easily 
understood: a trajectory of (56) for which KsmO* = 27rm with m integer, corresponds 
to a fixed point for 9 (because the angle is defined modulo 27r) and linear growth for J{t) 
(balhstic behavior). It can be shown that the stability region of these trajectories in phase 
space decreases as 1/K [52,30] and, for intermediate value of they play a important 
role in transport: particles close to these trajectories feel very long correlation times and 
perform very long jumps. The contribution of these trajectory, as a whole, gives the observed 
anomalous behavior. 
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Now, let us consider the cellular flow of Figure 11 as an example of sub-diffusive transport. 
We have seen that asymptotically (i.e. for t ^ ll/D) the transport is diffusive with effective 
diffusion coefficient which scales according to (64). For intermediate times Iq/vq ^ t <^ l^/D, 
when the boundary layer structure has set in, one expects anomalous sub-diffusive behavior 
as a consequence of fraction of particles which are trapped inside vortices [53]. A simple 
model for this problem is the comb model [31,28]: a random walk on a lattice with comb-like 
geometry. The base of the comb represents the boundary layer of size S around vortices and 
the teeth, of length /q, represent the inner area of the convective cells. For the analogy with 
the flow of Figure 11 the teeth are placed at the distance 6 ~ JDIq/vq (63). 



^0 

S 

FIG. 14. The comb model geometry 

A spot of random walker (dye) is placed, at time t = 0, at the base of the comb. In their 
walk on the x direction, the walkers can be trapped into the teeth (vortices) of dimension 
/q. For times Io/vq <^ t <^ l^/ D, the dye invades a distance of order {DtY/"^ along the teeth. 
The fraction F{t) of active dye on the base (i.e. on the separatrix) decreases with time as 

f W ~ ^ (69) 

and thus the effective dispersion along the base coordinate b is 

{b\t)) ^ F{t)Dt ^ 6{Dty/^ (70) 
In the physical space the corresponding displacement will be 

{x\t)) ^ {b\t))^j r-. k{PeDt)'/' (71) 
i.e. we obtain a sub-diffusive behavior with u = 1/4. 
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The above argument is correct only for the case of free-slip boundary conditions. In the 
physical case of no-slip boundaries, one obtains a different exponent = 1/3 [53]. The latter 
behavior has been indeed observed in experiments [17]. 

D. Beyond the diffusion coefficient 

Prom the above discussion it is now evident that diffusion, being an asymptotic behavior, 
needs large scale separation in order to be observed. In other words, diffusion arises only 
if the Lagrangian correlation time tl (60) is finite and the observation time is i ^ or, 
according to (62), if the dispersion is evaluated on scales much larger than Iq. 

On the other hand, there are many physical and engineering applications in which such 
a scale separation is not achievable. A typical example is the presence of boundaries which 
limit the scale of motion on scales L ~ Iq. In these cases, it is necessary to introduce 
non- asymptotic quantities in order to correctly describe dispersion. 

1 I ' -V ' ;, , ^ ' 1 

- 

-0.5 - ' ^ ^ - 

-1 ' ' ' 

-1 -0.5 0.5 1 

FIG. 15. Position of the 10000 tracers at a late stage of the point vortex advection dynamics. 

The position of the 4 vortices are in the center of the "hole" in which the tracers cannot enter. 

Before discussing the non- asymptotic statistics let us show, with an example, how it can 
be dangerous to apply the standard analysis in non-asymptotic situation. We consider the 
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motion of tracers advected by the two-dimensional flow generated by 4 point vortices in a 
disk. The evolution equation is given by (42) and (44) but now in (44), instead of Inr^, one 
has to consider the Green function G{rij) on the disk [39]. 

A set of 10000 tracers are initially placed in a very small cloud in the center of the disk. 
Because of the chaotic advection induced by the vortices, at large time we observe the tracers 
dispersed in all the disk (Figure 15). 

In the following we will consider relative dispersion, i.e. the mean size of a cluster of 
particles 

i?2(t) = (|x(t)-(x(t))|2) (72) 

Of course, for separation larger than the typical scale of the flow, Iq, the particles move 
independently and thus we expect again the asymptotic behavior 

R^{t) ~ 2Dt if R^itf^ » /o (73) 

For very small separation we expect, assuming that the Lagrangian motion is chaotic, 

R^{t) ~ i?2(o)e2^* if i?2^t)V2 ^ (74) 

where A is the Lagrangian Lyapunov exponent [23]. 

The computation of the standard dispersion for the tracers in the point vortex model 
is plotted in Figure 16. At very long time R^{t) reaches the saturation value due to the 
boundary. 

For intermediate times a power-law behavior with an anomalous exponent = 1.8 is 
clearly observable. Of course the anomalous behavior is spurious: after the discussion of the 
previous section, we do not see any reason for observing super-diffusion in the point vortex 
system. The apparent anomaly is simply due to the lack of scale separation and thus to the 
crossover from the exponential regime (74) to the saturation value. 
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FIG. 16. R^{t) computed for the four point vortex system. The dashed Hne is the power law 

To partially avoid this kind of problem, it has been recently introduced a new indicator 
based on fixed scale analysis [5]. The idea is very simple and it is based on exit time 
statistics. Given a set of thresholds (5„ = Sor'^, one measures the exit time Tj((5„) it takes 
for the separation Ri{t) to grow from 5n to Sn+i- The factor r may be any value > 1, but it 
should be not too large in order to have a good separation between the scales of motion. 

Performing the exit time experiment over particle pairs, from the average doubling 
time {T{5)) = l/NJ2iTi{S), one defines the Finite Size Lyapunov Exponent (FSLE) as 

In r 

- m) '''' 

which recovers the standard Lagrangian Lyapunov exponent in the limit of very small sep- 
arations A = lim^^o ^(^)- 

The finite size diffusion coefficient D{d) is defined, within this framework, as 

D{S) = 6^X{S) (76) 

For standard diffusion D{6) approaches the diffusion coefficient D (see (73)) in the limit of 
very large separations {S ^ Iq). This result stems from the scaling of the doubling times 
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{T{6)) ~ 5^ for normal diffusion. 

Thus, according to (73)- (74), the asymptotic behaviors of the FSLE are 



m 



(77) 



A if 5 < Zo 

if 5 > /o 

In presence of boundary at scales L ~ Iq, the second regime is not observable. For 
separation very close to to the saturation value Smax — L one expects the following behavior 
to hold for a broad class of systems [5] : 



X{6) oc 



(78) 



Let us now come back to the point vortex example of Figure 15. The FSLE for this 
problem is plotted in Figure 17. 
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FIG. 17. FSLE \{S) for the tracers advected by the point vortex system. The horizontal line 
represent the Lagrangian Lyapunov exponent (A ~ 0.14). The dashed curve is the saturation 
regime with Sjnax — 0.76. 

With the finite scale analysis one clearly see that only two regime survive: exponential 
at small scales (chaotic advection) and saturation at large scale. The apparent anomalous 
regime of Figure 16 is a spurious effect induced by taking the average at fixed time. 
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The finite scale method can be easily applied to the analysis of experimental data [11]. 
An example is the study of Lagrangian dispersion in a experimental convective cell. The 
cell is a rectangular tank filled with water and heated by a linear heat source placed on the 
bottom. The heater generates a vertical plume which induces a general two-dimensional 
circulation of two counter-rotating vortices. For high values of the Rayleigh number (i.e. 
heater temperature) the flow is not stationary and the plume oscillates periodically. In 

these conr^'^' — ^ — — 4,,™„ ^ j„ +^ +i, — f plume as a 

consequei 




FIG. 18. An example of trajectories obtained by PTV technique in the convective cell at 
Ra = 2.39 X 10^. The vertical thermal plume is clearly observable. The dark circle on the bottom 
represents the heat source. 

The study of Lagrangian dispersion has been done by means of the FSLE [11]. In 
Figure 19 we plot the result for Ra = 2.39 x 10*^. Again, because there is no scale separation 
between the Eulerian characteristic scale Iq (vortex size) and the basin scale L we cannot 
expect diffusion behavior. Indeed, the FSLE analysis reveals the chaotic regime X{S) = A at 
small scales and the saturation regime (78) at larger scale. 
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FIG. 19. FSLE X{S) computed for the convective cell at different initial separations (different 
symbols). The straight line is the Lyapunov exponent and the dashed curve represents the satura- 
tion regime. 

The finite scale tool has been successfully applied to many other numerical and experi- 
mental situations, from the dispersion in fully developed turbulence, to the analysis of tracer 
motion in ocean and atmosphere [33,13,38], to engineering laboratory experiments. It will 
be probably became a standard tool in the analysis of Lagrangian dispersion. 
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